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Abstract 

The thermo-magnetic corrections to the quark-gluon vertex in the presence of a weak magnetic field are calculated in the frame of 
the Hard Thermal Loop approximation. The vertex satisfies a QED-like Ward identity with the quark self-energy calculated within 
the same approximation. It turns out that only the longitudinal vertex components get modified. The calculation provides a first 
principles result for the quark anomalous magnetic moment at high temperature in a weak magnetic field. The effective thermo- 
magnetic quark-gluon coupling shows a decreasing behavior as function of the field strength. This result supports the observation 
that the behavior of the effective quark-gluon coupling in the presence of a magnetic field is an important ingredient in order to 
understand the inverse magnetic catalysis phenomenon recently observed in the lattice QCD simulations. 
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1> Introduction - 

Recent lattice QCD results, with 2-1-1 quark flavors, in¬ 
dicate that the transition temperature measured from the 
behavior of the chiral condensate and susceptibility, as 
well as from other thermodynamic observables such as 
longitudinal and transverse pressure, magnetization and 
energy and entrmy densities, decreases with increasing 
magnetic held Uslsl]- This result was called magnetic 
anticatalysis. 
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Recently, we have shown, in a natural way, how the 

decrease of the coupling constant with increasing held 
strength can be obtained within a perturbative calcula¬ 
tion beyond the mean held approximation. Notice that 
such a behavior cannot be achived in the mean held ap¬ 
proximation 10] ■ This was done both for the abelian 
Higgs lit] model as well as for the linear sigma model 
Ihl] where charged fields are subject to the effect of a 
constant magnetic field. This behavior introduces a de¬ 
pendence of the boson masses on the magnetic held in¬ 
ducing then a decreasing behavior of the critical tem¬ 
perature for chiral symmetry breaking/restoration. Go¬ 
ing now into QCD, in order to to establish if a similar 
behavior takes place, a first step would be to determine 
the finite temperature and magnetic held dependence of 
the coupling constant. 

The discussion of the thermal behavior of systems in¬ 
volving massless bosons, such as gluons, in the pres¬ 
ence of magnetic fields is quite subtle. Unless a careful 
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treatment is implemented, severe infrared divergences 
associated to the effective dimensional reduction of the 
momentum integrals will appear. These divergencs are 
associated to the separation of the energy levels into 
transverse and longitudinal directions (with respect to 
the magnetic field direction). The former are given in 
terms of discrete Landau levels. Thus, the longitudinal 
mode alone no longer can tame the divergence of the 
Bose-Einstein distribution. In this context, recently it 
has been shown that it is possible to find the appropriate 
condensation conditions by accounting for the plasma 
screening effects 13]. Here we show that a simple pre¬ 
scription where the fermion mass acts as the infrared 
regulator allows to obtain the leading behavior of the 
QCD coupling for weak magnetic fields at high tem¬ 
perature, that is, in the Hard Thermal Loop (HTL) ap¬ 
proximation. As a first step, we compute the thermo- 
magnetic corrections to the quark-gluon vertex in the 
weak field approximation. We use this calculation to 
compute the thermo-magnetic dependence of the QCD 
coupling. To include the magnetic field effects we use 
Schwinger’s proper time method. We should point out 
that the weak field approximation means that one con¬ 
siders the field strength to be smaller than the square of 
the temperature but does not imply a hierarchy with re¬ 
spect to other scales in the problem such as the fermion 
mass. 


2. Charged fermion propagator in a medium 


The presence of a constant magnetic field breaks 
Lorentz invariance and leads to a charged fermion prop¬ 
agator which is function of the separate transverse and 
longitudinal momentum components (with respect to 
the field direction). Considering the case of a magnetic 
field pointing along the z direction, namely B = Bz, the 
vector potential, in the so called symmetric gauge, is 

Af,{x) = ^(0,-X2,xuO). ( 1 ) 

The fermion propagator in coordinate space cannot 
longer be written as a simple Fourier transform of a mo¬ 
mentum propagator but instead it is written as 

S(x,x')^<b(x,x') J £^e-'>(^-^'^S(k), (2) 


where 


<l>(x, x') - exp 



+ -F^y(^ - x'r 


( 3 ) 

2 


is called the phase factor and q is the absolute value 
of the fermion’s charge, in units of the electron charge. 
S (k) is given by 


S(k) 


- -‘f 


ds 




cos{qBs) 

X I [cos(qBs) - 1 - 7172 sin(g'Bs)] (m -i- ^y) 


\ 

cos(qBs))’ 


( 4 ) 


where m is the quark mass and we use the definitions for 
the parallel and perpendicular components of the scalar 
product of two vectors a^ and given by 


(a ■ b)\\ - aobo - a^bj, 

(a-b)± - aibi+a 2 b 2 . (5) 


P2 P2 



Figure 1. Feynman diagrams contributing to the thermo-magnetic 
dependence of the quark-gluon vertex. Diagram (a) corresponds to 
a QED-like contribution whereas diagram (b) corresponds to a pure 
QCD contribution. 


Figure [T] shows the Feynman diagrams contributing to 
the quark-gluon vertex. Diagram (a) corresponds to 
a QED-like contribution whereas diagram (b) corre¬ 
sponds to a pure QCD contribution. The computation 
of these diagrams requires using the fermion propaga¬ 
tor given by Eq. (|2]i, which involves the phase factor in 
Eq. (O. It turns out, however, that this phase does not 
contribute. We refer the reader to the original article 
where a proof of this fact is presented. Therefore, 
for the computation of diagrams (a) and (b) in Fig. [1] 
we can just work with the momentum representation of 
the fermion propagators since the phase factors do not 
contribute. The situation would have been nontrivial in 
case the computation had required a three fermion prop¬ 
agator closed loop M- 

3. QCD vertex at finite temperature with a weak 
magnetic field 

To compute the leading magnetic field dependence 
of the vertex at high temperature, we work in the 
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weak field limit of the momentum representation of the 
fermion propagator We work in Euclidean space 
which is suited for calculations at hnite temperature in 
the imaginary-time formalism of finite temperature field 
theory. The fermion propagator up to 0{qB) is written 
as 


S{K)^ 


m — 

+ irP- 


- iy 172 


(K^ + irP-y 


iqB). 


( 6 ) 


Using the propagator in Eq. (Hll and extracting a factor 
gta common to the bare and purely thermal contribu¬ 
tions to the vertex, the magnetic field dependent part of 
diagram (a) in Eig.[T]is expressed as 


= -igHCf-CJ2)(qB)T^ J ^ 

X Tv [rir2^iir/i^A(^2 - k) 

+ ¥^7^l7l72K\A(Pl - 7v 

X A{K)A{P2 - K)A(Pi - K), (7) 


where, in the spirit of the HTL approximation, we have 
ignored terms proportional to m in the numerator and 

1 


A{K) = 


A(K) = 


+ irP 

1 




l+k^ 


( 8 ) 


with a)n - (2n + 1 )7tT and tj„ - 2nnT the fermion and 
boson Matsubara frequencies. Cf, Ca are the factors 
corresponding to the fundamental and adjoint represen¬ 
tations of the S U{N) Casimir operators 


two channels are already accounted for in Eq. O since 
the magnetic held insertion on each quark internal line 
is thereby included at the order we are considering. 
These expressions can be simplihed by noting that in 
the HTL approximation the external momenta Pi and 
P 2 are small and can be taken of the same order. We 
have carried out the calculation of the Matsubara sums 
but we do not have space here to present all the details. 
Adding both contributions to the vertex corrections we 
hnd a quite compact expression 

= 2iy5g^CF{qB)G^{Pi,P2), ( 11 ) 

where 


G,,(PuP2) = 


X 


2 ry 

vj (271^ 

|(/: ■ bWu^ - (K ■ uWb^ 


+ [(K ■ b)^i- (K ■ 

X A^(K)A(Pi - K)A(P2 - K). (12) 


It is no surprise that the thermo-magnetic correction to 
the quark-gluon vertex is proportional to 75 since the 
magnetic held is odd under parity conjugation. It is 
important to stress that our vertex satishes a QED-like 
Ward identity, i.e. 


(El - P 2 ) ■ 6 r(Pi,P 2 ) = E(Ei) - E(E 2 ), (13) 


Ca = N, (9) 

respectively. Hereafter, capital letters are used to re¬ 
fer to four-momenta in Euclidean space with compo¬ 
nents Kfj - {kn, k) = (-Cl), k), with w either a Matsubara 
fermion or boson frequency. 

Eollowing the same procedure, the magnetic held de¬ 
pendent part of diagram (b) in Eig.[T]is expressed as 



X \-¥y\y2¥\\yii + 2rvrir2^iiry^/j 

- 7/^7172^11^] 

X 'A(KfA(Pi - K)A(P2 - K). (10) 

The explicit factor 2 on the right-hand side of Eq. (fTOl i 
accounts for the two possible fermion channels. These 


where E(E) denotes the self energy correction for the 
quark in the presence of an external magnetic held. This 
can be computed from the diagram in Eig. |2] and the 
corresponding expression in the HTL-approximation is 
given by 

9 1 r cTk 

^(P)-2iysg^CAqB)T^ J ^ 

A(P - K)A\K) xI^(K ■b)^-(K ■ M)|i|. (14) 

The technical details requeried for computing the fre¬ 
quency sums cannot be presented here due to lack of 
space. The reader is refered once again to |@]. 

This remarkable result shows that even in the presence 
of the magnetic held and provided the temperature is the 
largest of the energy scales, the thermo-magnetic cor¬ 
rection to the quark-gluon vertex is gauge invariant. 
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Figure 2. Feynman diagram for the the quark self-energy. The inter¬ 
nal quark line represents the quark propagator in the presence of the 
magnetic field in the weak field limit. 


4. Thermo-magnetic QCD coupling 

The calculation of G^{P\,P 2 ) demands not only to deal 
with Matsubara sums but also to discuss tensor expres¬ 
sions of the type 


JaiiPuPl) 


I 


dO. kaki 
47r {P,.k){P2-k)' 


(15) 


For the sake of simplicity we have chosen a particular 
configuration where the external momenta P\ and P 2 
make a relative angle 6*12 = n. It turns out that only the 
longitudinal components of our tensor, i.e of the ther- 
momagnetic vertex, are modified. If we rearrange the 
gamma matrices to introduce the spin operator in the z 
direction S 3 , we find for the explicit longitudinal com¬ 
ponents 


^f||(po) = 


' 2 ' 


Ag^CpM^iT, m, qB)yi^2, 


(16) 


where Jn - (yo, 0 , 0 , - 73 ) and where we have defined 
the function M^(T, m, qB) as 


M^{T, m, qB) 


qB 

\67jP- 


ln( 2 )-^- 


(17) 


That the vertex correction in Eq. (fThl l is proportional to 
the third component of the spin operator is of course nat¬ 
ural since the first order magnetic correction to the ver¬ 
tex is in turn proportional to the spin interaction with the 
magnetic field, which we chose to point along the third 
spatial direction. The correction, thus, corresponds to 
the quark anomalous magnetic moment at high temper¬ 
ature in a weak magnetic field. Note also that Eq. (O 
depends on the scales po ™d m. po is the typical en¬ 
ergy of a quark in the medium and therefore the simplest 
choice for this scale is to take it as the temperature. The 
quark mass represents the infrared scale and it is there¬ 
fore natural to take it as the thermal quark mass. We 
thus set 


Po ^ T 


The pure thermal correction to the quark-gluon vertex is 
well known lEl] to be given by 


6rf^^(PuP2)--mj f 


dQ k^k 

{p,-k){P2-k)' 


(19) 


In order to extract the effective modification to the cou¬ 
pling constant in one of the longitudinal directions, 
let us look at the contribution proportional to yo from 
Eq. (fT9l l. Eor the same working momenta configuration 
we find 

2 

lYi 

<‘'"™(fo) = ~yo. (20) 

Po 

Using Eqs. (fTSl l and (l20l i. the effective thermo-magnetic 
modification to the quark-gluon coupling extracted from 
the effective longitudinal vertex, is given by 


^eff = 8 


iri f / 8 \ 

I - + \-^yCFM\T,mf,qB) 


, ( 21 ) 


Eigure [3] shows the behavior of ggff normalized to 
^therm> where 


^ therm — 8 



( 22 ) 


for Qs - g^l4-7T - 0.2,0.3 as a function of the scaled 
variable b - qB/T^. Plotted as a function of b 
and for our choice of momentum scales, the function 
^eff/^therm temperature independent. 



Figure 3. The effective thermo-magnetic coupling normalized to 
the purely thermal coupling as a function of the field strength 

scaled by the squared of the temperature, for = 0.2,0.3. Note that 
ggff decreases down to about 15% — 25% for the largest strength of 
the magnetic field within the weak field limit with respect to the purely 
thermal correction and that the decrease is faster for larger 
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Note that the effective thermo-magnetic coupling ggp- 
is a decreasing function of the magnetic field. The 
decrease becomes more significant for larger as and 
for the considered values of as it becomes about 15% 
-25% smaller than the purely thermal correction for 
qB~T^ ~ 1 . 

5. Summary and Conclusions 

In this work we have computed the thermo-magnetic 
corrections to the quark-gluon vertex for a weak mag¬ 
netic field and in the HTL approximation. It turns out 
that this vertex satisfies a QED-like Ward identity with 
the quark self-energy. This results hints to the gauge- 
independence of the calculation. The thermo-magnetic 
correction is proportional to the spin component in the 
direction of the magnetic field, affecting only the longi¬ 
tudinal components of the quark-gluon vertex. It thus 
corresponds to the quark anomalous magnetic moment 
at high temperature in a weak magnetic field. 

From the quark-gluon vertex we extracted the behavior 
of the magnetic field dependence of the QCD coupling. 
The coupling decreases as the field strength increases. 
For analytic simplicity, the explicit calculation has been 
performed for a momentum configuration that accounts 
for the conditions prevailing in a quark-gluon medium 
at high temperature, namely namely back-to-back slow 
moving quarks whose energy is of the order of the tem¬ 
perature and with the infrared scale of the order of the 
thermal particle’s mass. Notice that the chosen values 
for these scales provide a good representation of the 
coupling constant’s strength within the plasma condi¬ 
tions. 

We stress that the weak field approximation means that 
the field strength is smaller than the square of the tem¬ 
perature but does not require a hierarchy with respect to 
other scales in the problem such as the fermion thermal 
mass. Note that for situations where the temperature is 
not the largest of the energy scales, as for the largest 
magnetic field strengths achieved in peripheral heavy- 
ion collisions, or for compact astrophysical objects, the 
present calculation does not provide conclusions about 
the behavior of the coupling constant nor for its effect 
on the critical temperature. Nevertheless, it is impor¬ 
tant to note that lattice QCD simulations for the critical 
temperature are performed for magnetic field intensities 
starting from qB - 0 (see e.g. Fig. 9 in Ref. lEl]). The 
results of these simulations show that the critical tem¬ 
perature is a decreasing function of the field intensity 
all the way to qB = I GeV^. Also, one should note that 
it is likely that the effects on the quark-gluon plasma in 
peripheral heavy-ion collisions take place when the field 


intensity is already smaller than the square of the tem¬ 
perature, since the field strength is a rapidly decreasing 
function of time. It is thus important to explore in QCD 
the regime where qB is small compared to as a nec¬ 
essary bridge to further studies on the behavior of the 
coupling constant for larger field strengths. 

The result supports the idea iSi that the decreasing 
of the coupling constant may be one of the important 
ingredients to understand the inverse magnetic catalysis 
obtained in lattice QCD. 
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